The Kondo lattice model from strong-coupling 

viewpoint 

Valeri N. Kotov* and Peter J. Hirschfeld 

Department of Physics, University of Florida, Gainesville, FL 32611-8440, USA 
Abstract 

We study the magnetic excitation spectrum of the two-dimensional (2D) square- 
lattice S=l/2 Kondo lattice model at finite (hole) doping, by representing the 
Hamiltonian in terms of a set of local excitations. The location of the paramagnetic- 
antiferromagnetic phase boundary at T=0 is determined. 
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The Kondo lattice model in its simplest form describes itinerant (c) electrons 
hopping on a lattice, and interacting with localized (d) orbitals via Kondo 
exchange J (we assume an antiferromagnetic sign J > 0): 

H = -tY: {cUc^, + h.c.) + JE S - -S^ (1) 

i,j,a i 



There are two major issues that are of interest from both fundamental and 
experimental (heavy fermion materials) perspectives: (1.) Description of the 
competition between the Kondo effect and the RKKY interaction, which gen- 
erally leads to a quantum phase transition from a paramagnetic to a long-range 
ordered phase, and (2.) The possible breakdown of the Kondo effect (screen- 
ing) near such a phase boundary, especially in antiferromagnetic metals [1]. 
We concentrate mostly on the first aspect, and consider the model Eq.(l) on 
a square lattice at finite hole doping 5 (5 = is half-filling), which makes the 
system conducting. 

We rewrite Eq.(l) in the local (one site) basis, where four bosonic modes 
(a singlet and a triplet) can be excited: V2s^v) = e,p44|0), V2tl\v) = 
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ap^esacldj^lO) , (a = x,y,z; (J,p,S =t,i), as well as two fermionic modes: 

fi,a\'^) = "^ctIO); /l.c-l"^) = cjc|(ij.|0). Here \v) is a reference vacuum with no 
excitations present, and |0) corresponds to absence of physical electrons. The 
described set of states can be used to generalize the bond-operator representa- 
tion, often used in localized spin systems (with only s and t modes present) [2]. 
Similar ideas have also been successfully applied to Kondo insulators {S — 0) 
[3]. The representation is exact provided the on-site kinematic constraint is 
implemented: sKs + fl^fu,cr + ^l^a = 1) necessary to restrict the Hilbert space 
of the new operators to the physical sub-space. It is physically clear that for 
doping 5 — Q and strong Kondo exchange J/t ':$> 1 the ground state is a 
singlet (thus paramagnet), meaning that the s excitations condense: (s) ^ 0. 
Let us note that the term "strong-coupling" (J/t ^ 1) in this work has the 
opposite meaning to the one used in Ref.[l]. The Hamiltonian in the local 
representation becomes, by taking the singlet state s^v) as the ground state, 

H = H2 + Hint- 

H,^ Y: £;.(k)/t,,k/.,.,k + JEtk-tk + E^(k)/iV,k/2U 

k,(T,f=l,2 k k,<T 

Hint = E [tt//] ^ [t//] scattering vertices, (3) 

where i?jy(k), -B(k) oc t are easily calculated, and Hint represents the non-linear 
interactions (with coefficients also proportional to the hopping t) between 
the bosonic and fermionic modes. The challenge now becomes to treat the 
kinematic constraint (which amounts to an infinite local repulsion between the 
t and fi, modes), as well as the interactions present in Hint, as accurately as 
possible in order to have reliable results for the excitation spectrum. We have 
followed the philosophy outlined in Ref . [4] for spin insulators to implement the 
hard-core repulsion, while the various terms in Hint were treated to lowest non- 
trivial (one- loop) order (see also references [5] and [6] , discussing the situation 
at finite doping for the t — J model) . The extension of the approach from zero 
to finite doping is a non-trivial matter, and details will be provided elsewhere 
[7]. It is crucial in this treatment that the S=l magnetic excitations (the 
spectrum of the t operators) form a dilute Bose gas. Departure from the dilute, 
low-density limit would also mean a possible change in the ground state and 
failure of our approximation scheme since the latter is based on classification 
of diagrams in powers of the excitation density. In addition, doping is also 
assumed to be a small parameter in the problem. 

Let us summarize our main results [7]: 

• The magnetic (t) excitations are found to be gapped everywhere in a certain 

region of parameter space, see Fig.l, meaning that the system is paramagnetic. 
The gap at the antiferromagnetic ordering wave vector Qaf = (^r, vr) vanishes 
on a critical line, signaling a transition to a phase with long-range order (fi- 
nite staggered magnetization, proportional to the condensate {tk=QJ^p)). The 
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parameter F, displayed on the axis in Fig.l (which is introduced for numerical 
reasons), is found to vary very weakly with doping. In the simplest of our 
approximation schemes it has the constant value Ffti 0.72 [7]. This produces a 
critical point at zero doping (J/t)c ~ 1-39 which is quite close to the numerical 
result [8] (however no numerical results are available for 5 7^ 0). Let us also 
mention that the magnons are damped at any finite doping, which can readily 
be seen from the width of their spectral function. 

• We would like to emphasize that our calculational scheme is self-consistent 
in the sense that the dilute Bose gas description is maintained throughout the 
phase diagram. This means that the singlet ground state, the starting point 
of the calculation, is stable (i.e. (s) ^ 0), leading to the usual Fermi hquid de- 
scription of the critical properties [9] . Departure from the Fermi liquid picture 
would occur if (s) ^ at the quantum critical point (QCP); one can then 
expect various exotic effects such as different critical exponents, which has 
been in fact suggested to occur in two dimensions [1] . The fact that we do not 
find perturbative indications for breakdown of the Kondo effect at the QCP 
could be quite possibly due to our approximation scheme, which can detect 
the existence of a quantum critical line (where [the magnetic gap]<^ J), but is 
not capable of predicting universal properties in the quantum critical regime. 
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Fig. 1. Phase diagram of the Kondo lattice model. 
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